Long-range sound-mediated dark soliton interactions in trapped atomic condensates 
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A long-range soliton interaction is discussed whereby two or more dark solitons interact in an 
inhomogeneous atomic condensate, modifying their respective dynamics via the exchange of sound 
waves without ever coming into direct contact. An idealized double well geometry is shown to 
yield perfect energy transfer and complete periodic identity reversal of the two solitons. Two 
experimentally relevant geometries are analyzed which should enable the observation of this long- 
range interaction. 

PACS numbers: 67.85.-d, 03.75.Lm, 05.45.Yv. 



I. INTRODUCTION 

Solitons arc an important phenomenon present in 
many areas of nonlinear physics [l[ arising from the bal- 
ance between dispersion and nonlinearity. Dark solitons 
are localized waves propagating within a background 
medium, and are characterized by a dip of the density 
and a phase slip 0] ■ They have recently become a topic 
of intense research in weakly-interacting atomic Bosc- 
Einstein condensates Q, being routinely generated in a 
number of experiments, both in a controlled manner [34^1 
and as a result of dynamical processes (Iol - fl4l |. 

Recent experiments have been able to generate appro- 
priate low-temperature conditions in quasi-lD geometries 
enabling the observation of one or more (undamped) dark 
soliton oscillations in a harmonically-confined condensate 
[ll-dl as predicted at the mean field level [L5l4l7| . In 
general, the presence of the axial confinement in these 
experiments breaks the integrability of the system, ren- 
dering the soliton unstable to sound emission along its 
axis of propagation, inducing it to decay [TH fl8l. Il9j]. 
Earlier work by us has revealed the crucial role of the 
harmonicity of the underlying trap in stabilizing the 
propagating soliton against decay at low temperatures 
via continuous cycles of soliton-sound interactions [H, 
|20| , when all other decay mechanisms are suppressed for 
the timescales of interest [UrM H3411] . 

Here we show that soliton-sound interactions also play 
a key role in the dynamics of two or more solitons in an 
inhomogeneous atomic condensate, and we identify opti- 
mal realistic experimental conditions for observing such 
an effect. In particular, we analyse the motion of solitons 
oscillating within either a single harmonically-confined 
condensate or different spatially-separated sub-regions of 
a condensate and conclude that, within mean field theory, 
dark solitons can interact via the emission/absorption of 
sound waves over a long range, without ever approaching 
very close to each other. We refer to this interaction 
mechanism as sound-mediated dark soliton interactions, 



and suggest an experiment that could probe such an 
effect. 

Although this effect bears close analogies to earlier 
studies in various physical systems and configurations 
(see below) I33l - l36| | , this study is distinct in that (i) 
the integrability of the system is lifted by the longitudinal 
inhomogeneous confinement, as in current cold atoms ex- 
periments, and that (ii) at sufficiently low temperatures 
such exchange need not necessarily lead to the decay 
of the solitons involved, but it can instead lead to a 
periodic exchange of energy between them. We wish to 
note here that direct collisions of two dark solitons were 
previously studied both in a ho mog eneous [liH l34j and a 
harmonically confined medium [351 ] . and where found to 
result at most in a phase shift of the emerging solitons. 
Long-range dark soliton interactions have also been 
discussed in [33| , in the context of an axially homogeneous 
condensate, in which solitons interact via superradiance 
due to the quantised nature of the background fluctua- 
tions. Finally, the proposed effect is reminiscent of long- 
range bright soliton interactions observed in non-local 
nonlinear optical media over distances much larger that 
the soliton size 1361. 
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FIG. 1: (color online) Motion of two non-identical dark 
solitons in a harmonically trapped condensate. Left: Density 
(black) and potential (dashed blue). Right: Position of 
left (solid, grey) and right soliton (solid, black) vs. time. 
(Initial parameters: zl = — 3.0Z Z , v s = 0.3c and Zr = 0.0Z Z , 
vf = 0.1c, (jl = 25hto z ). 
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II. DEMONSTRATION OF EFFECT 
A. Two solitons within a single harmonic trap 

Previous work 0, [H| found that two identical, grey 
solitons displaced equally from the bottom of the har- 
monic trap and travelling towards each other, arc re- 
flected or pass through each other depending on their ki- 
netic energies, indicating that soliton-sound interactions 
have, at most, a minor effect on their mutual dynamics. 
Although of very small magnitude, such modulations 
are actually present even in that context; however, this 
effect does become significant in the case of non-identical 
solitons initially located at asymmetric locations within 
a single harmonic trap, as shown in Fig. 1. 

As well-known, if there is just one soliton propagating 
in a harmonic trap, it perturbs the inhomogeneous 
background motion, setting up a dipole mode, but it 
also periodically reabsorbs the sound emitted by its 
own motion [181 ] : the resulting oscillation exhibits a 
beating effect due to the different oscillation frequencies 
of the soliton (= w x / V2) [H-ii, HI [H] and the induced 
background dipolar motion of the condensate (= u> z ). 
In the case of two non-identical solitons with different 
initial speeds, the solitons are found to absorb sound and 
dissipate at different rates, as each soliton re-interacts 
both with the sound it has emitted and with sound 
emitted by the other soliton. This process, combined 
with periodic soliton collisions, leads to a much more 
pronounced modulation of the soliton oscillations than 
encountered in the case of a single soliton in a harmonic 
trap (whose corresponding amplitude modulation would 
typically be a few percent of its initial amplitude (37|V 

This effect is shown clearly in Fig. [1] (Right): in 
particular, the fast soliton, initially imprinted away 
from the trap centre, gains energy, thereby exhibiting 
oscillations of decreasing amplitude during the initial 
time period, t < 70W" 1 . This energy becomes available 
from the sound emission originating from the motion of 
the slower soliton, which in turn exhibits a small increase 
in its oscillation amplitude; at later times, the change in 
the relative phase of the two soliton oscillations leads to a 
reversal of the direction of the energy flow between them. 
Clearly, however, their mutual energy exchange does not 
lead to a net decay within a single harmonic trap. 

Our analysis is based on the one-dimensional (ID) 
Gross-Pitaevskii Equation 

<-(-^ +v « + 'i*i'-")*' (1) 

where g parametrizes the effective ID interatomic in- 
teraction via g = (4nh' 2 a/m)/27rl 2 L , where a is the 
3D s-wave scattering length; here l± = \Jfij (mw_i_) is 
the harmonic oscillator length in the tightly-confining 
transverse direction, of effective harmonic frequency u>±, 
m is the atomic mass and il is the chemical potential. For 
the trap parameters of a recent experiment used here Q , 



this correspond to an atom number of the order of 1000 
atoms. 



B. Two solitons in a spatially separated idealized 
double-harmonic trap 

To isolate the two effects, direct soliton collisions and 
soliton-sound interaction, in a manner that could also 
be achieved experimentally, we now move to a modified 
geometry which enables us to restrict the motion of the 
two solitons in two spatially separated traps, while still 
allowing sound exchange between them. Due to the 
known sensitivity of soliton-sound interaction on any 
trap anharmonicity [2(|, we initially choose an idealized 
geometry V(z) = mui^{z ± zq) 2 /2, with two minima 
at ±zo, such that the region z<0 provides a 'harmonic 
environment' for each soliton individually; this physical 
separation enables us to probe soliton-sound interactions 
over a long spatial range, and for long times. 

The amount of sound transferred between the two 
sides of the trap is set by the ratio of Vq/li, where 
Vq = tow^Zq/2 is the height of the barrier between the 
neighbouring harmonic traps. Throughout this work, 
we fix this ratio to Vo/fx ~ 0.9, which allows sound 
transfer between the two wells while keeping the solitons 
individually confined in different wells. 

To monitor this idealized energy exchange mechanism, 
we consider for simplicity a stationary black soliton in the 
centre of the right well (zr = 7l z ), and a moving, 0.3c 
(grey) soliton in the centre of the left well (zl = —7l z ), 
and track their subsequent coupled evolution. As the left 
soliton oscillates, it emits sound; while part of it remains 
confined to the left well, thus re- interacting with the same 
soliton, the remaining part moves across the barrier into 
the other well (since Vq < li). The sound entering the 
right well disturbs the initially stationary soliton present 
there, and causes it to move, itself radiating sound, 
part of which also becomes transferred into the left well. 
This process of emission, rcabsorption and mutual sound 
exchange between the two solitons continues and results 
in a periodic reversal of their characteristics. This effect 
is clearly shown in Fig. [5] (bottom). One can think of 
one soliton ac ting as a parametric external driver for the 
other soliton [381 ]. with the frequency of the driver set 
dynamically by soliton-sound interaction. 

Note that when the soliton in one well is moving the 
fastest, its amplitude of oscillation is at a maximum, 
and the soliton in the other well is stationary (its 
depth is maximum), and vice versa. This transfer of 
energy is more apparent when plotting the position 
of the soliton in the left well against the position of 
the soliton in the right well at different times of their 
coupled evolution, as shown in Fig. [2] (top) . This figure 
corresponds to a temporal duration of approximately one 
soliton oscillation in the trap, i.e. r = 2V2TTUJ- 1 ; the 
bottom figure shows the individual soliton oscillations, 
highlighting the times around which the top images were 
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FIG. 2: (color online) Top images: Snapshots of soliton 
position in the left trap, zl, against that in the right trap, 
zr, at various times. Bottom: Oscillations of two solitons 
(«s = 0.3c, v, — 0) in an idealized double harmonic trap 
( with Vo/n = 0.98, zl = — 7l z , zr = 71 z ), showing the 
complete interchange of the identities of the two solitons. Blue 
lines indicate the growing oscillation envelope of the initially 
0.3c speed soliton placed in the left well when there is no 
soliton in the other well. 



taken. Note that after a few more cycles the maximum 
amplitude of oscillation will gradually start increasing, 
resulting in an imperfect energy transfer; eventually both 
solitons will decay (even within pure mean field theory), 
due to the anharmonicity in the region of the barrier. 
Increasing Vb/V will result in observation of the identity 
reversal over a much longer period of time, conversely, 
reducing the ratio will enhance it. However, reducing it 
too much will result in a more pronounced anharmonic 
region and so the solitons will be prone to acceleration 
and dissipation will occur in a shorter timescale. 

This behavior should also be compared to that of a 
moving soliton in the left well, when the part of the 
condensate trapped in the right well does not contain a 
soliton. Unlike the case of a dark soliton in a 'closed' 
harmonic trap geometry, the presence of a region for 
the emitted sound to escape leads to a slow decay; the 
envelope of such anti-damped oscillations are shown by 
the blue lines in Fig. 2 (bottom). The comparison to 
this scenario demonstrates clearly that the presence of 
the soliton in the right well acts so as to stabilize the left 
soliton, enabling it to survive/oscillate for much longer 



III. EXPERIMENTAL CONSIDERATIONS: 
DOUBLE/MULTIPLE WELLS 

We seek to identify a directly observable effect that 
could be measured in experiments. Despite the appealing 
feature of the 'identity-reversal' of the two solitons in 
the idealized cut-off harmonic trap discussed previously, 
such an effect is unlikely to be directly observed in an 
experiment; this is due to the sensitivity of solitons to 
the locally strongly anharmonic region near the barrier 
separating the two regions in a realistic experimental 
double-well geometry. As an alternative, and to make 
direct links with experiments, we therefore consider the 
simpler scenario of looking at changes in the motion of a 
soliton in one well, due to the presence, or absence, of a 
soliton in the neighbouring well. 

Thus, we now discuss experimental geometries in 
which solitons have been successfully generated, namely 
Gaussian traps and optical lattice geometries. For 
simplicity we restrict this study to identical solitons. 



A. Gaussian induced, double-well potential 

We first consider a double well geometry of the form 
V(z) = mojlz 2 /2 + Voe~ z / 2d ; such a geometry can 
be formed by adding a repulsive Gaussian dimple to a 
harmonic trap, e.g. by the addition of a blue detuned 
laser [39j . with Vq and d characterizing the height and 
width of the Gaussian. As before, two solitons are placed 
at the peak densities of the two wells spanned by the 
condensate, with their respective profiles shown in Fig. [3] 
(Left). 

The wells have been chosen such that the effective 
harmonic frequency of each well is of the order of the 
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FIG. 3: (Color Online) Left: Density profile (black and grey 
dashed) and trapping potential (lower dashed blue) of an 
atomic condensate in a double-well potential, with either a 
single 0.1c soliton in the right well (zr = 28. 5£, black lines), 
or also a second identical soliton located in the symmetric 
position of the left well (zl = 28. 5£, green dashed). Right: 
Position of the soliton in the right well when it is the only 
soliton present in the trap (black) and when there is another 
soliton in the left well (grey). (Parameters: Vo/fi = 0.97, 
d = 14.6£.) 
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longitudinal frequency u> z used earlier; however, within 
each well, the pronounced anharmonicity which arises 
from the lack of symmetry leads to the decay of the 
soliton. Hereafter we scale distances to the healing length 
£ (£ * W4.9 for the previous parameters). 

Given the rapid decay of a soliton in such an asymmet- 
ric trap, we instead investigate the motion of a soliton in 
the right well depending on whether there is a soliton in 
the well to the left or not. We find that the addition of 
a soliton in the left well stabilizes the other soliton to 
a large extent, as shown in Fig. [3] (Right); as a result, 
instead of quickly decaying, it continues oscillating with 
roughly the same amplitude (i.e. same energy/depth) 
for a timescale of the order of 800£/c, before eventually 
decaying. We also note a slight phase difference between 
the two cases. 



B. Optical double- well potential 

To increase the soliton lifetimes we next probe a regime 
in which each well of the trap is itself symmetric, but also 
not harmonic. Intersecting laser beams at an angle can 
be used to generate an optical lattice or an array of traps 
of variable depth and periodicity. This configuration is 
modelled by V(z) = Vocos 2 (2nz/d), with the lattice 
parameters Vb and d characterizing the lattice's depth 
and spacing respectively. Again, solitons of equal speed 
are initially placed at the peak condensate densities. 

As before, we initially consider a periodic array of 
traps of relatively large periodicity [40], in which the 
condensate however only spans two wells, as shown in 
Fig. 0] (top) . We investigate the evolution of an initially 
propagating soliton (v = 0.1c) in one well, when the 
condensate in the other well does not contain a soliton 
(black line); this is contrasted to the case when there 
is an identical soliton in the other well (grey line). It 
is apparent that the presence of the other soliton leads 
to a slower decay. This indicates that the soliton in the 
right well absorbs the sound created by the motion of the 
soliton through the inhomogeneous background of the left 
well. The more symmetric nature of the trap results in 
the soliton experiencing a slower dissipation than in the 
Gaussian, double- well potential of Fig. EH 



C. Array of optical traps 
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FIG. 4: (Color Online) Left: Density profiles for a condensate 
spanning (a) two, or (c) six sites of an array of optical traps. 
In both cases, the right site contains a v 3 — 0.1c soliton 
(black). In addition, in (a) the left site may or may not 
contain an identical soliton (green dashed), whereas in (d) 
the leftmost site also contains an identical soliton, while the 
condensate in the intermediate sites may or may not contain 
solitons. Right: Temporal evolution of the soliton in the 
right well showing clearly that the presence of a soliton in 
neighbouring wells has a stabilizing effect: (b) Comparison of 
case when there is no soliton in the left well (black) to when 
this site contains a soliton (grey), (d) Comparison of the 
dynamics of the soliton in the rightmost trap, when there are 
no intermediate solitons (black) and when the intermediate 
sites contain identical solitons (grey). (Parameters: Vo/n = 
0.87 and d = 170.9£). 



the partial escape of the emitted sound to the other wells, 
while the sound from the soliton in the leftmost well 
arrives too late to stabilize it. This should be contrasted 
to the case when each intermediate site contains solitons 
of the same depth. In that case, the solitons in the 
neighbouring wells clearly have the net effect of partially 
stabilizing the end soliton against decay over a prolonged 
period, as evident by the grey lines in Fig. 4(d). 

Faster solitons were found to decay on a more rapid 
timescale, as the time required for the oscillating soliton 
to probe the anharmonic region of the trap decreases 
rapidly with increasing soliton speed. However, this 
effect should be observable in carefully controlled exper- 
iments, and it does not require identical soliton speeds. 



Finally, we wish to explicitly demonstrate the role of 
dark solitons as both absorbers and emitters of sound [l8| , 
by investigating the effect of adding further solitons to 
a condensate spanning a number of lattice sites greater 
than two. 

We consider a condensate spanning 6 lattice sites, and 
place two identical slowly-moving solitons in the wells at 
each end of the condensate, as shown in Fig. 4 (bottom). 
The soliton in the right well is found to decay (black lines 
in Fig. 4(d)) on a timescale of ~ 25 oscillations, due to 



IV. CONCLUSIONS 

In conclusion, we have demonstrated that two or more 
solitons propagating in an inhomogeneous medium can 
interact and exchange energy without ever coming into 
direct contact, or colliding with each other; the physical 
mechanism of this exchange is based on the mutual 
exchange of sound energy over distances much larger than 
the soliton size. One could also interpret this process as 
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follows: one soliton perturbs the background condensate, 
setting up a dipole mode in it; this dipole mode, perturbs 
the condensate region where a second soliton is confined, 
thereby forcing the other (propagating) soliton to see a 
temporally-varying background condensate; this leads to 
a change in the interaction of the second soliton with the 
background, which can periodically lead to absorption of 
energy from the background, rather than emission. In an 
idealized geometry, this interaction leads to a complete 
energy transfer which results in periodic identity reversal 
of the two solitons. 

Under current experimental conditions this identity 
reversal may not be observable, as more oscillations are 
required than have been observed up until now 
and one would additionally need to significantly minimize 
the region where the potential becomes anharmonic. 
However, experiments should be able to detect that the 
presence of a second soliton in a condensate confined in 
an appropriate double well geometry has a stabilizing 



effect. The experiment we have in mind requires an 
experimental set up which can reproduce the same 
initial conditions to good accuracy, so that differences 
in individual soliton trajectories in one well after a 
number of 10-15 oscillations arising from the presence of 
another soliton located in a spatially separated well can 
be observed (e.g. after expansion imaging). Generalizing 
this to an array of traps (which could be thought of 
as a soliton gas [4l|) the addition of more solitons was 
found to lead to a slower dissipation of the outer solitons. 
Extensions of this work could include the interplay 
between soliton-sound interaction and Josephson effects 
[42| , a deeper investigation of the role of the intermediate 
states in a multi-site configuration in terms of resonant 
tunnelling [43l |44| , and the role of quantum and thermal 
fluctuations. 
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